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The  second  order  Doubly  Asymptotic  Approximation  (DAA2) 
i im  t Tn nl  ilmeli>]iinl  hr  flosrs  and  r<  lip|x  (I'fdm  been  implemented 
in  the  SWEEPS  STAGS  codi-f#tfo  provide  4  steady -state  vibration 
analysis  capability  (or  submerged  structures  This  paper  extends 
the  previous  work  irf  PtRnnis.  fltmri  nrt  P»npps  |li?l  1r  deter¬ 
mining  the  convergence  behavior  of  DAA2-predtctsd  wet  surface 
variables  that  are  relevant  to  steady -stale  submerged  structural 
vibrations  The  convergence  rates  of  *^ffuid  boundary  eigenvalues"^ 
and  convergence  behaviors  of  modal  acoustic  impedance  functions 
are  calculated  for  increasingly  dense  fluid  element  grids  on  a 
spherical  surface  Within  mode-dependent  frequency  ranges  In 
which  DAA2  is  valid,  dense  grids  on  the  order  of  five  elements  per 
naif  wave  are  needed  for  five  percent  accuracy  In  modal  acoustic 
impedance  This  IS  a  consequence  of  the  low  order  (constant 
source  strength)  interpolation  used  for  establishing  fluid  variables 
There  are  mode-dependent  frequency  ranges  for  which  DAA2 
acoustic  impedance  predictions  are  inaccurate,  even  for  very  dense 
fluid  element  grids  These  results  provide  guidance  for  engineers 
developing  SWEEPS  STAOS  idealizations  for  vibration  analysis  of 
complex  submerged  structures^ 

NOMENCLATURE 

lOverdots  denote  differentiation  with  respect  to  time,  e  g  .  p  is 
first  time  derivative  of  pressure) 


1  -  total  number  of  fluid  elements  in  spherical  surface  mesh 


acoustic  wave  number 


structural  suffneas  matrix 


m  -  order  of  asymmetric  spherical  surface  harmonic 
l ml  -  superscript  denoting  mode-derived  DA A2 
m.  -  nth-mode  wet  surface  flu.d  mass 


wet  surface  fluid  masa  matrix 


structural  mass  matrix 


order  of  axUymmetrlc  spherical  surface  harmonic 
wet  surface  pressure 

wet  aurfaoe  pressure  vector  associated  with  fluid  boun¬ 
dary  eigenmode 

incident  wav*  wet  mir face  pressure  vector 


nth  mode  wet  surface  element  area 


wet  surface-radiated  pressure  vector 


A.  -  wet  surface  element  area  matrix 


spherics!  surface  radius 


fluid  sound  speed 


Laplace  transform  variable 


ici  *  superscript  denoting  curvature-augmented  DAA2 


wet  surface  normal  velocity 


Cs  -  structural  damping  matrix 

f,  •  mechanical  input  force  vector 

g  -  mode-derived  DAA2  empirical  curvature  parameter 
G  -  structure  coordinate  -  fluid  coordinate  iransfomauon 


wet  surface  normal  velocity  vector  associated  with  fluid 
boundary  eigen  mode 

incident  wave  wet  surface  normal  velocity  vector 
nth  mode  wet  surface  normal  velocity 
wet  surface  normal  velocity  vector 


g  -  structure!  displacement  qgtsr 

Zj,  -  nth-mode  (specific)  »flomitc  fluid  nwniaw 

k  •  diagonal  matrix  of  wet  aurfaea  maan  local  curvature# 

K„  -  nthmoda  wat  surface  mean  curvature 

1  -  diagonal  matrix  of  fluid  boundary  tifsnvaluee 

1„  »  nth  moda  fluid  boundary  alganvalue 

e  ■  fluid  danalty 

Q,  ■  wat  surfaos  fluld-frequancy  matrix 

u  -  radian  frequency 

un  -  nth-mode  wet  surface  fluid-frequency  parameter 


INTRODUCTION 

the  steady-jtate  time- harmonic  interaction  between  submerged 
elastic  structures  and  exterior  acoustic  fluid  Is  an  important  Naval 
structural  mechanics  and  acoustics  problem.  The  finite  element 
method  (FEM)  la  the  only  viable  means  of  modeling  the  vibratory 
reepor.ee  of  harmonically -excited,  topographically  oomplax  struc¬ 
ture  oompoeed  of  beams,  plates,  and  eball*  The  FEM,  aa  a  means 
of  approximating  the  maaa,  stiff neaa.  and  damping  diatrlbutfona  of 
arbitrary  structure*.  la  daaorlfaad  in  many  taxta;  [4]  la  a  well- 
known  example.  Various  numerical  echemee  are  available  for 
"exact"  modeling  of  Linear  aoouatlc  fluid  in  tlma- harmonic  pro¬ 
blems;  two  important  approaches  are  fluid  finite  elements  (8.6,7] 
and  HMmhqlts  aurfaoe  integral  equations  [BJJ.lQl.  Both  of  thaas 
fluid  models  have  been  coupled  to  structural  finite  aliment  oodee 
(11-15)  such  that  vibration  and  aound  radiation  problama  of 
coupled  fluid-structure  dynamical  ay  stems  can  be  analysed.  Fluid 
finite  element  and  Helmholtz  surface  integral  approaches  ars 
three-dimensional  and  two-dimensional  numerical  mathoda  of 
dlacretlzing  the  waee  equation.  Thaas  methods  are  capable  of  ex¬ 
actly  representing  Ould  loading  afreets  on  the  wat  aurfaoe  of  a 
submerged  vibrating  structure.  Convergence  toward  exact  solu¬ 
tion*  with  increasing  wat  surfaos  grid  density  is  posalbla  aa  long 
aa  propagating  acoustic  waves  in  the  fluid  domain  and  Quid  boun¬ 
dary  pressure  distributions  are  resolved  (Also.  In  otrtaln  Helm¬ 
holtz  integral  methods  the  problem  of  uniqueness  of  solution  has 
to  be  dealt  with). 

An  alternative  approximate  acoustic  fluid  model  is  the  Second 
order  Doubly  Asymptotic  Approximation  (DAAJ2).  originally 
developed  for  transient  fluid-structure  interaction  [16],  and  re 
oently  extended  to  uma-harmonlc  vibrations  of  submerged  struc¬ 
tures  (1).  DAA2  has  also  been  exercised  in  simple  coupled  vibra¬ 
tion  and  sound  radiation  problems  (17).  Mathematically .  DAAR  la  a 
seoond-order  differential  equation  that  approximates  the  exact 
relationship  hetwaan  fluid  pressure  end  velocity  on  a  two- 
dimensional  boundary  of  a  three  dimensional  lnflnita  fluid  do- 
main  The  boundary  equations  era,  affectively,  dimaniionally- 
roduoed  versions  of  three-dimensional  fluid  field  aquations 
Numerical  implementation  of  DAA2  follows  a  boundary  element 
method  (BEM.  also  discussed  in  (4)).  since  the  spatial  variations  or 
fluid  pressure  and  velocity  on  the  boundary  of  an  Infinite  exterior 
fluid  l wet  surface  of  a  submerged  structure)  ars  interpolated  by 
assumed  functions  local  to  elements  of  the  surface.  The  frequency 
domain  impedance  (ratio  of  pressure  to  velocity)  approximation  la 
governed  by  the  order  of  the  DAA  equation  used 

A  structural  dynamics  analyst  usually  invokes  many  approx¬ 
imations  whsn  devising  idealizations  of  complicated  structures 
The  inherent  assumptions  of  toe  mathematics  underlying  an 
analysis  computer  code  (e  g  .  structural  FEM  coupled  with  fluid 
BEM  in  the  preeent  ease)  comprise  an  approximation  lava)  to  be 
accepted  at  the  onset.  The  analyst  must  endeavor  to  minimus  the 
influence  of  discretization  errors  on  predicted  response  variables. 
More  often  than  not.  the  analyst  faces  the  task  of  predicting  the 


response  af  a  eubmargsl  structure  far  which  no  previous  modeling 
ekperlenoe  exists  In  such  caste,  information  on  the  ability  of  the 
chosen  numerical  method  to  represent  the  response  of  simple 
structures  demonstrating  known  behavior  la  useful  guidance  for 
developing  valid  idealisations  of  oompUoettd  structures 

Valid  FEM/BEM  Idealizations  of  vibrating  submerged  struc¬ 
tures  are  those  which  oorrectly  resolve  structural  bending  and 
membrane  travelling  waves  and  resonant  standing  waves  poorly 
oontrotlod  by  structural  damping-  Stmilarty.  the  aoouatlc  im¬ 
pedance  that  exterior  fluid  presents  to  the  structural  wat  surface 
must  be  resolved  in  the  frequency  domain  and  In  apace.  Mathe¬ 
matical  idealizations  have  hope  of  meeting  thaas  requirements  if 
they  possess  (1):  a  sufficient  number  of  structural  elements  par 
wavelength  of  important  banding  or  membrane  response  patterns, 
and  (6):  a  sufficient  number  of  fluid  boundary  elements  par 
wavelength  in  important  wat  aurfaoe  stationary  or  travailing 
response  patterns.  The  optimum  wat  aurfaoe  FEM/BKM  grid 
depends  on  the  Interpolation  orders  of  atructui  il  and  fluid 
elements  relative  to  the  response  pattern  complexity. 

This  paper  presents  benchmark  problem  results  that  should  be 
useful  to  user*  of  the  8WKKP8-8TAGB  cod*  (8.18).  which  utilize* 
coupled  structural  FEM  and  DAA2  baaed  fluid  BEM  for  prediction 
of  the  nonlinear  transient  or  linear  tlme-harmonlr  response  of 
submerged  shell  atructurea.  In  particular,  the  ability  of  the 
SWEEF8-STAQS  fluid  boundary  elements  and  the  DAA2  approxima¬ 
tion  to  resolve  the  modal  aoouatlc  impedance  functions  of  a 
spheric*!  surface  exposed  to  an  lnflnita  extant  of  external  acoustic 
fluid  la  examined. 

A  DAAL2  MODEL  FOR  ACOUSTIC  IMPEDANCE 

The  DAA2  equations  defining  an  approximate  relationship  be¬ 
tween  pressure  end  velocity  on  the  aurfaoe  of  e  submerged 
vibrating  structure  have  been  reported  elsewhere  (2|.  they  ere 
time-harmonic  specialization*  of  the  transient  equations  USi  The — 
coupling  eg  structural  stiffness,  structural  damping,  structural 
man.  and  complex-valued  acoustic  fluid  Impedance  Is  apparent  in 
the  DAAS  frequency  domain  matrix  equations  leqe.  Ill  end  |2>  of 
(21): 


s- 

Kft  ]  1  P.  > 

i«f) 

+  l"€.+  5. 

? .1  ’  ?A r 

-  oo 

Ef,  •  -w*Hf  ♦  aC  CJtoA,  +  »(  A,) 

«•“!*-  Q&rPi 

it  *  o°  <<“*¥1  -  MMfc>  Si 

An  approximate  model  of  the  pressure  end  velocity  relation 
ship  on  a  closed  surface  within  an  exterior  fluid  of  infinite  extent, 
shown  in  Figure  1.  can  be  obtained  by  suppressing  eii  structural 
metric**  end  specifying  zero  incident-wen*  pressure  and  velocity 
Under  these  conditions,  the  scattered  pressure  con  be  interpreted 
as  radiated  pressure  and  equation  ill  reduces  to 

1  -w*M,  +  pc  (JoA,  +  0,A,l|  pr  +  lee  <J<**U|CT  +  <U*0,M,CTI|  A  =  0  i2) 

The  following  transformation  converts  "structural  displace 
ment  vector  g  to  wet  surface  velocity  vector  u. 


u,  »  ->u>Gtx 

If  g„  is  imagined  aa  describing  modal  standing  wave  motions 
of  the  fluid  "boundary"  in  the  absence  of  external  driving  forces 


3 


Uhj.Vifk  £f»  ViS-tWVi'Vv'V  v.%  CVVf. .."t.;  UV»  AW-—; 


i.vtWitvi'.i-:- 


bwAw- 


V*  W  a.“-  '*'*  V  •  w  •— »  V>  V 


TWO-WMWWWWAL  IURFACI 


INFINITE 
EXTERIOR  aUD 


Figure  1  -  Two-Dimensional  Fluid  Surface  In  Thrae-Dlmenslonal 
Exterior  Fluid  Domain 

fs.  then  (2).  when  combined  with  (3),  reduce  to  a  frequency- 
domain  relationship  between  fluid  boundary  pressures  and 
velocities  required  to  sustain  these  modes: 

l-cvaM*  -*■  yc  I ja>Ar  +  QfA{)|  p  a  [pc(JwQfMf  -  (naMf)|  u  (4) 


modal  lmpedanoe  is  a  function  of  one  fluid  boundary  mode,  eaoh 
of  whloh  oorreeponde  to  the  low  frequency  limit  of  the  exaot  "ac- 
omlon  to  inertia  per  unit  area”  of  an  acouetlo  mode  of  a  epherloal 
surface. 

The  general  expressions  for  DAA™1  and  DAA'|’  modal  acoustic 

Impedances,  equations  (A16)  and  (A17)  of  the  Appendix,  reduce  to 
simpler  forms  when  the  normalizations  e  “  c  ”  r  -  unity  are  used 

for  a  spherical  surface  geometry.  The  DAA™1  empirical  curvature 
parameter  g  Is  taken  to  be  unity  for  the  spherical  surface,  and  the 
DA  A1  g1  curvature  function  Qn  reduces  to  1-An.  With  these  particu¬ 
lars,  equations  (A16)  and  ( A  1*7)  become: 

(oil,,)4  -  «*>*„)*  +  1  (6) 

C _ 

(<uAn)4  -  (wA,/  +  1 

_ 

(wAn)4  -  «»*„>*  +  8<V 

,  (7) 

twAn>4  -  <OiAn>*  +  g(An) 
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An  equivalent  expression  In  the  time  domain  Is: 


where  gUn)  -  (A*  +  1)  +  2Jn((o»An)‘  -!)■ 


Mf  p  +  pcAf  p  +  fcQrAr  p  a  pc  |M,U  +  QrM(  u)  (6) 

Equations  (4)  and  (5)  are  the  DAA2  equations  defining  the 
fluid  pressure  p  on  the  surface  of  a  fluid  boundary  vibrating 
freely  with  velocity  distribution  u  at  frequency  w.  This  Interpreta¬ 
tion  of  the  conditions  under  which  the  general  DAA2  fluid- 
structure  expressions  (1)  were  reduced  to  either  (4)  or  (8)  Implies 
that  DAA2  modal  acoustic  Impedance  functions  can  be  defined  for 
’very  fluid  boundary  "mode.''  provided  the  modes  are  uncoupled. 
This  Is  possible  only  for  a  limited  number  of  fluid  surface 
geometries;  the  Infinitely  long  cylindrical  surface  and  the  spherical 
surface  are  two  examples.  The  spherical  surface  Is  considered 
here. 


The  "exact"  low  frequency  limits  of  a  spherical  surface's  fluid 
boundary  eigenvalues.  An  -  l/(n+l),  can  be  used  In  equations  161 
and  (7)  to  define  "exact"  DAA2  modal  acoustic  impedance  func¬ 
tions  Zn.  Approximate  An,  as  predicted  for  various  8WEEPS-8TAGS 
fluid  element  grids,  can  be  substituted  In  the  same  expressions  to 
obtain  associated  approximate  Zn.  The  convergence  of  Zn  and  A„ 
toward  values  based  on  exact  fluid  boundary  eigenvalues  (the  best 
results  possible  with  DAA2)  with  Increasing  fluid  element  grid 
density  is  now  examined.  Frequency  ranges  of  DAA2  method 
validity  are  also  determined  by  comparison  to  exact  modal  series 
solutions  [21). 

SWEEPS-STAGS  MODAL  IMPEDANCE  RESOLUTION  8TUDV 


SPECIALIZATION  TO  SPHERICAL  SURFACE 

The  acoustic  Impedance  functions  Implied  by  the  DAA2  dif¬ 
ferential  equations  relating  wet  surface  pressure  and  velocity  are 

derived  for  the  two  DAA2  forms  (mode-derived  DAA™'  and 

curvature-augmented  DAA™)  in  the  Appendix.  Equations  (4)  and  (5! 

above  are  used  as  a  basis  for  this  development.  When  these  two 
functions  are  specialized  to  uncoupled  acoustic  mode  situations, 
the  resulting  complex-valued  DAA2  acoustic  Impedance  poly¬ 
nomials  are  functions  of  "fluid  boundary  mode”  eigenvalues  and 
radlated-wave  curvature  parameters.  The  notion  of  fluid  boundary 

modes  and  the  meanings  of  DAA™'  and  nAA™  are  described  more 
completely  by  fleers  [16.191  and  Fellppa  [20).  In  DAA™'  the  curv¬ 
ature  parameter  is  a  mode-independent  empirical  constant,  while 
DAA'?,'  uses  a  matrix  of  local  wet  surface  curvatures.  As  shown  In 

the  Appendix,  the  accuracy  of  DAA2-predlcted  wet  surface  modal 
acoustic  impedance  depends  on  the  accuracy  of  computed  fluid 
boundary  eigenvalues  Fluid  olement  grid  densities  needed  for  a 
given  acceptable  acoustic  impedance  error  level  In  frequency 
ranges  of  DAA2  method  validity  can  be  readily  defined  for  the 
spherical  wet  surface  problem  since  the  modal  impedances  are  un¬ 
coupled  As  shown  In  the  Appendix,  eaoh  DAA2-based  complex 


Equations  (6)  and  (7)  show  that  DAA2  polynomial  approxima¬ 
tions  of  Impedance  functions  Zn  depend  only  on  frequency  w  and 
fluid  boundary  eigenvalues  An.  Thus,  only  the  FLUMAS  (Fluid 
Mass)  module  of  the  SWEEPS-STAGS  code,  which  calculates  fluid 
mass  matrix  and  solves  the  following  fluid  boundary  eigen- 
problem,  needed  to  be  run  In  this  study: 

M,u  =  AAfu  i8> 

The  fluid  boundary  eigenvalues  and  eigenvectors  of  a  spherical 
surface  of  unit  radius.  Immersed  In  a  fictional  acoustic  fluid  of 
unit  density  and  sound  speed,  were  determined  for  a  series  of 
uniform  SWEEPS-STAGS  fluid  element  meshes.  The  mesh  family 
considered  Is  shown  in  Figure  2:  the  four-noded  fluid  elements  of 
the  FLUMAS  module  of  SWEEPS-STAGS  were  used  In  all  cases. 

As  shown  In  DeRuntz  and  Geers  [3).  a  three-dimensional 
spherical  surface  possesses  repeating  eigenvalues  of  many  nonax- 
isymmetrlc  wave  number  orders  (m)  for  every  value  of  axlsym- 
metric  wave  number  order  (nl  of  n  =  2  and  above.  Hence.  FLUMAS 
produced  as  many  repeated  eigenvalues  for  n  >  2  as  were  allowed 
by  the  symmetry  plane  constraints  specified  In  the  quarter  sphere 
idealizations  of  Figure  2.  These  eigenvalues  occurred  In  the  se 
quence  (1  ln+11).  which  is  a  reduction  of  equation  8.31  of  dunger 
and  Felt  |22)  for  unit  spherical  surface  radius  and  unit  fluid 
density. 
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Figure  2  -  Family  of  Fluid  Element  Meshes  on  a  Spherical  Surface 
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Figure  3  •  Convergence  of  Spherical  Surface  Fluid  Boundary 
Eigenvalues  with  Increasing  Fluid  Element  Density 


In  (3),  convergence  rates  of  fluid  boundary  eigenvalues  were 
presented  for  n=0  and  n=l  modes  of  a  spherical  surface.  The  pres¬ 
ent  work  extends  these  results  to  n  >  2.  FLUMAS-computed  fluid 
boundary  eigenvalue  convergence  rates  for  the  grids  of  Figure  2 
are  presented  In  Figure  3  These  computations  exhibit  a  conver¬ 
gence  rate  of  roughly  (1/1),  where  I  Is  the  number  of  elements  In 
the  mesh  These  predictions  agree,  qualitatively,  with  those  obtain¬ 
ed  earlier  by  DeRuntz  and  Oeers  [3],  The  predicted  differences  bet¬ 
ween  theoretically  Identical  eigenvalues  of  like  n-order  but  differ¬ 
ing  m-order  become  greater  as  n  Increases  and/or  grid  density 
decreases  This  Is  Indicated  by  the  error  ranges  drawn  In  Figure  3. 


co/ln+l) 


Figure  4  -  Erroneous  Mode-Derived  DAA2  Prediction  of  Acoustic 
Impedance  of  N  -  0  Mode 
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Figure  6a  -  Comparison  of  Sweeps-Stags  Modal  Acoustic  Impedance 
Predictions  to  Exact  Solution  for  Spherical  Surface 

The  DAA2  spherical  surface  Impedance  functions  of  equations 
10)  and  (7)  were  calculated  for  each  axlsymmetrlc  mode  0  <  n  <  4, 
and  for  each  mesh  of  Figure  2  In  the  normalized  frequency  range 
0.1  <  o>/(n+l)  <  10.0.  It  Is  known  that  acoustic  Impedance  follows 
a  transition  between  high-  and  low-frequency  asymptotic  limits  in 
this  frequency  range.  The  results  of  the  Impedance  calculations 
are  shown  In  Figures  4  and  6. 

The  mode-derived  forms  of  DAA2  modal  acoustic  Impedances 
given  by  equation  (6)  were  found  to  be  Inadequate.  Although  con¬ 
vergence  with  Increasing  grid  density  occurred,  the  converged 
functions  were  totally  erroneous.  Figure  4  demonstrates  this  fund¬ 
amental  Inadequacy  of  DAA1™  by  comparison  to  exact  fluid  im¬ 
pedance  functions  of  Fellppa  and  Geers  [21). 

In  certain  frequency  ranges,  the  DAA  ^  polynomials  of  equa¬ 
tions  (7),  with  exact  fluid  boundary  eigenvalues  substituted,  did 
not  accurately  represent  the  exact  modal  series  acoustic  impedance 
solutions  [21).  This  error  reflects  DAA2  frequency  domain  approxi¬ 
mations  and  is  Independent  of  fluid  element  grid  density.  On  the 
other  hand,  acoustic  Impedances  based  on  approximate  fluid  boun¬ 
dary  eigenvalues  did. converge  to  the  exact  modal  series  solution  In 
frequency  ranges  where  the  exact  DAA2  solutions  are  fair  repre¬ 
sentations  of  exact  modal  series  acoustic  Impedance. 

The  different  convergence  behaviors  for  n=0  through  n=4  are 
Illustrated  in  Figure  5.  The  modal  series  n=0  impedances  are  ac¬ 
curately  predicted  by  DAA1^'  In  the  entire  frequency  range  0.1  < 
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Figure  Sc  Comparison  of  Sweeps  Stags  Modal  Aoousuc  Impedence 
Predictions  to  Exact  Solution  Tor  Spherical  Surface 
loonunuedi 
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Figure  Se  Comparison  of  Sweepe  Stags  Modal  Aoousuc  Impedance 
Predictions  to  Exact  Solution  for  Spherical  Surface 
i  continued) 


w  in*  1 1  o  10  0  even  with  a  very  oosrse  grid  (Figure  5a>  The 
DAA  '  n-1  impedancea  converged  to  the  exact  fluid  boundary 
eigenvalue  soluuon  with  increasing  grid  density.  For  n-1.  the  ex 
act  DAA  £  soluUons  are.  like  for  n»u.  essenuaiiy  the  exsct 
modal  sene*  curves  (Figure  Elbl  However  for  n-2.  3 
and  4.  the  exact  DAA1'1  solutions  are  no;  ulus:  to  the  modi!  scries 
solutions  in  the  entire  frequency  range  of  interest  Exact  modal 
series  exact  DAA  '  and  approximate  DAA  n-2  modal  acoustic 

impedances  are  compared  in  Figure  5c  for  several  fluid  element 
grids  For  n-2.  the  frequency  domain  errors  dominate  the  acoustic 
reactance  dm  Z,,l  in  the  nondimensional  frequency  range  0  4  < 
,„,in-*li  <4  0.  and  dominate  the  acoustic  resistance  (Be  Zn)  for  w 
<  1  4  This  means  that  DAA2  is  inherently  mcapatle  of  properly 
representing  the  "added  mesa"  effect  m  n-2  vibrations  of  a 
spherical  surface,  as  governed  by  Im  Z(i.  in  the  dimensional  fre¬ 
quency  range  12  <u/  <  12  0.  and  also  will  not  accurstely  simulate 
radiation  damping  for  u,  <  4  2  Figures  5d  and  5e  show  similar 


limiting  ranges  for  n*3  and  4  The  frequency  ranges  of  validity  for 
DAA  prediction  of  n-2.  3.  and  4  modal  acoustic  impedances  are 
summarised  In  Tsbie  l 

The  next  consideration  is  the  reduction  of  spatial  discretiza¬ 
tion  error  in  the  applicability  ranges  listed  in  Table  1.  The 
predicted  magnitudes  of  impedance  terms  for  the  several  grids 

considered  here  are  compared  to  modal  series  and  exact  r/AA 

solutions  in  Tables  2  and  3  Another  limiting  facto*  illustrated  .n 
these  tablet  It  the  inability  of  coarser  fluid  element  grids  to 
resolvs  fluid  boundary  eigenvalues  and  recognizable  eigenvectors 
The  ratios  of  predicted  to  txvct  DAA2  fluid  boundary  eigenvalues, 
(low  frequency  asymptotic  values  of  accession  to  inertia  per  unit 
area)  are  gi-’en  for  various  grid  danslltet  in  Table  4  FLU  MAS  did 
predict  multiple  eigenvalues  for  the  nonaxisymmelric  modal  pat¬ 
terns  allowed  by  the  two  symmetry  planes  in  the  quarter  sphere 
wet  surface  idealizations  In  a  similar  calculation  by  DeRuntz  i  18). 
a  384  element  grid  was  unable  to  predict  multiple  eigenvalues  for 
the  iionaxisymmetnc  m-11  modal  order  For  this  waveform  the 
erroneous  differences  betwoen  predicted  degenerate  eigenvalues 
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TABLE  1  •  FREQUENCY  RANGES  OF  VALIDITY  FOR  CURVATURE- 
AUGMENTED  DAA2 


TtOS  OF  Zn  COMPONENTS:  (EXACT  DAA'j'/ EXACT  MODAL  SERIES 


^■xuiKri'AiKEUiKruiKasKixail 


TABLE  3  •  CONVERGENCE  OF  DAAj,  /  MODAL  ACOU8TIC 
RESISTANCE  FOR  VARIOU8  SPHERICAL  8URFACE  GRIDS 


c.uaixi?idaK 


1  /EXACT  MODAL  SERIKSl 


0.971  0.894 

0.999  1.000 


ERROR  IS  OVER  FIVE  PERCENT  IN  RANGES: 


<u  <  1.4 
Cl)  <  1.4 
W  <  1.6 


lm  Zn 

(REACTANCE) 


R«  Zn 

(RESISTANCE) 


0  986 
1 .002 
1.009 
1.012 
1.012 


(*)  :  MESH  INADEQUATE  TO  DEFINE  FLUID  BOUNDARY  EIGENVALUES 
(*)  :  OUTSIDE  RANGE  OF  DAA^  VALIDITY. 


1.004 

1  004 

1.022 

1  005 

1  028 

1  005 

TABLE  2  •  CONVERGENCE  OF  DAA^  7  MODAL  ACOUSTIC  REACTANCE 
FOR  VARIOUS  SPHERICAL  SURFACE  GRIDS 


RATIOS  OF  lm  2n  IDAA  j  7  EXACT  MODAL  SERIES) 


ml 

m 

ifffi 

SBE 

HRS 

P'P 

csa 

SEkHc 

■  *714* 

LU 

MTTTW 

LLUl 

KLil 

KXI1S 

■mi 

wtuzm 

mW'y.m 

1X3 

0  857 

0  960 

0  981 

0  988 

0993 

0  978 

0  994 

1  004 

1.017 

01 

0  848 

0.968 

0  980 

0984 

0  993 

0  979 

0  999 

1009 

1.023 

0  4 

0  854 

0  982 

0  943 

0  949 

0  993 

0  958 

0  979 

0.989 

1004 

04 

0  885 

0  975 

0  989 

0  993 

0  196 

0  0 

0  942 

0*93 

09*8 

om 

0  990 

0  7 

1  014 

1  013 

1  007 

1  004 

1.003 

0  B 

1  088 

1  029 

1  014 

1.009 

1.005 

0  9 

1  142 

1  040 

1  018 

1  011 

1.007 

1*1 

1  0 

1  179 

1  046 

1  021 

1  013 

1  008 

1  2 

1  203 

1  048 

1.021 

1  013 

1  008 

1.4 

1  192 

1  042 

1  019 

1.010 

1  007 

I  6 

1  189 

1  038 

1.018 

1  009 

1  006 

11 

1  148 

1  03C 

1  014 

1  007 

1  005 

2  0 

1  125 

"" 

1  012 

1  003 

1  004 

_ 

■mu 

H 

Ki 

knv^h 

0  1 

' 

0.984 

1005 

0  2 

0.996 

0  3 

<*l 

u  ■ 

1.008 

0  4 

Bnl 

0  997 

0  5 

■ 

0  945 

...  MESH  INADEQUATE  TO  DEFINE  FLUID  BOUNDARY  EIGENVALUES 
l.l  OUTSIDE  RANGE  OF  QBA<C|  VALIDITY 


were  of  the  same  order  as  the  spacing  between  eigenvalues  of  adja¬ 
cent  axlsymmetric  n-orders.  All  of  this  data  Implies  minimum  re¬ 
quired  grid  densities  for  acoustic  resistance  and  reactance  predic¬ 
tion  In  frequency  ranges  of  DAA1'1  validity,  as  limited  by  ln 
resolution. 

Table  5  summarizes  the  convergence  behavior  of  n=l,  2,  3  and 
4  fluid  boundary  eigenvalues,  acoustic  reactances,  and  acoustic 
resistances  for  several  different  SWEEPS-STAGS  fluid  element 
grids.  The  number  of  eloments  per  half-wave  of  the  modal  pattern 
required  to  resolve  fluid  boundary  eigenvalues  and  modal  acoustic 
impedance  functions  accurate  to  five  percent  depends  somewhat  on 
the  fluid  eigenmode  shape.  The  n=l  mode  seems  to  require  special 
consideration  It  Is  fair  to  conclude,  however,  that  a  minimum  of  5 
elements  per  half-wave  Is  needed  for  five  percent  accuracy  In  fluid 
boundary  eigenvalue  and  resistive  fluid  Impedance  prediction  In 
the  range  u.  >  1.4  for  n=2  and  3  and  for  <u  5.  1.6  for  n=4.  Similar 
requirements  seem  to  hold  for  minimizing  spatial  discretization  er¬ 
rors  til  reactive  impedance  ("addod  mass"  offoat)  predictions  at 
relatively  low  frequencies.  These  tentative  mesh  requirements  In 
SWEEP  STAGS  submerged  structure  vibration  analysis  may  be 
substantially  altered  by  Inclusion  of  structural  Impedance  and/or 
structural  damping  Also,  convergence  of  far-fleld  radiation  will 
probably  not  bo  is  sensitive  to  grid  density  as  the  near-flold 
variables  considered  In  this  work 


TABLE  4  -  CONVERGENCE  OF  DAX^'FLUID  BOUNDARY 
EIGENVALUES  FOR  INCREASING  SPHERICAL  SURFACE 
GRID  DENSITY 


RATIOS  OF  DAA(§'  A„:  (APPROX. /EXACT) 

m 

E29 

1=24 

1=96 

1=210 

1=600 

0 

VO 

»3 

w 

M3 

996 

541 

1 

.50 

.874 

.965 

.984 

990 

994 

2 

.3333 

.798-  88 

.946 

976 

MS 

Mi 

3 

.26 

820-  .88 

.924.932 

".964” 

9B6 

MB 

4 

.20 

.876  1.015 

.910  .925 

955 

.975 

980 

THEORETICALLY  DEGENERATE  MULTIPLE  EIGENVALUES  OF  LIKE  n  ORDER.  DIFFERING 
m- ORDER.  NOT  RESOLVED  FOR  n.  1  COMBINATIONS  BELOW  SOLID  LINE 


TABLE  6  -  REQUIRED  SWEEPS-STAGS  FLUID  ELEMENT  MESH 
DENSITIES  FOR  ACCURATE  PREDICTION  OF  MODAL 
ACOUSTIC  IMPEDANCE 


ELEMENTS  PER  MODAL  HALF  WAVE 

jjgjt 

n  =  1 

rv  2 

n  3 

n  4 

2 

4 

6 

8 

HI 

HALFWAVES 

HALFWAVES 

HALFWAVES 

HALFWAVES 

mm 

8 

4 

2 

1.33 

1 

96 

16 

_ JL  J 

4 

2.67 

2 

216 

24 

12 

6 

3 

384 

32 

16 

8 

5  33 

4 

600 

40 

20 

10 

6  67 

5 

MESH  DENSITIES  BELOW  THESE  LINES  RESOLVE  THESE  VARIABLES  TO  WITHIN  FIVE  PERCENT  0E 
EXACT  VA.UES 


.  FLUID  BOUNDARY  EIGENVALUES  A„ 

- MODAL  ACOUSTIC  REACTANCE  lm  Z„ 

- MODAL  ACOUSTIC  RESISTANCE  Re  Z„ 

CONCLUSIONS 

This  work  provides  some  initial  guidelines  on  required 
SWEEPS-STAGS  fluid  grid  density  needed  for  accurate  representa¬ 
tion  of  acoustic  Impedance  on  a  surface  exposed  to  exterior  fluid. 
These  guidelines  should  be  useful  In  developing  analytical  Idealiza¬ 
tions  of  complicated  submerged  structures.  The  required  fluid  grid 
density  for  five  percent  error  in  acoustic  Impedance  Is  five 
elements  per  half-wave  In  frequency  ranges  whore  DAA2  Is  valid. 
This  requirement  Is  rather  demanding  and  certainly  roflocta  the 
low  Interpolation  order  of  fluid  variables  In  the  SWEEPS-STAGS 
fluid  elements.  Further  studies  should  be  undertaken  to  determine 
whether  the  total  wet  surface  Input  Impedance  of  a  submerged 
structure  of  Interest  Is  governed  by  structural  Impedance  (mass, 
stiffness,  and  structural  damping!  or  by  acoustic  Impedance.  It  Is 
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likely  that  the  *u irfeo*  grid  requirement*  for  uiuttqi  of  strue- 
tunl  vibration  prediction  would  vary  eooordirw  to  Uw  topographic 
details  of  the  structure,  the  distribution  and  phasing  of  the  Input 
loads,  and  lha  stiff  nee*  of  Local  rations  of  the  structure  On  the 
other  hand,  accurals  prediction  of  fsr-fltld  radlalsd  sound 
pras«urs  m  the  fluid  remote  from  the  structure  would  probably 
not  be  as  aensltiv*  to  surface  (rid  density 
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APPENDIX 

MODAL  IMPEDANCE  FUNCTIONS  IMPLIED  BY  D*A2 
DIFFERENTIAL  EQUATIONS 

The  DAA2  differemisl  equation  infatmg  **»  vrt*  surfs re  rstc 
city  vector  u  to  lha  surface  pressure  vector  p  is.  for  coupled  modal 
responses 

M(p  +  pcA,p  +  ecB,A,p  -  pc  (M,u  +  @,M(U)  (All 

The  matrix  g,  is  on  explicit  function  of  wet  surface  curvature, 
and  is  an  implicit  function  of  the  wet  surface's  "fluid  boundary 
modes"  which  are  the  eigenfunction*  of  the  zero-frequency  hydro¬ 
dynamic  alganproblem: 

Mfu  •  !A,u  1*2) 

In  "mode  derived"  DAA2  (DAA’g 9f  »  defined  by  the  follow¬ 
ing.  with  the  curvature  parameter  "g"  restricted  to  the  given 
range 

S’7'  “  gecAfM,'1.  0«  g<  1  0  lA3i 

In  "curvstinw-sijgmented"  DAA2  fSAA'^'l.  O,  is  defined  in 
terms  C  l  matrix  of  local  wet  surface  curvature*  k 


-  fcAfM,-'  -  ck 


In  DAA  '™\  the  parsmeur  g  i*.  heuristicelly  s  measure  of  the 
average  curvature  of  acoustic  waves  emitted  by  the  radiaung  wet 
eurfaee  The  lower  bound  value  g«0  reduces  second  order  DAA'™’ 

to  first  order  DAA1.  in  which  the  limiting  high  frequency  fluid  im¬ 
pedance  asymptote  la  pc.  the  impedance  exhibited  by  plane  wave* 
The  other  extreme.  g“  1  0.  ha*  produced  reasonable  results  in  s 
USA-STAGS  aphencal  shell  scattering  problem  (18).  while  g-0  5 
was  found  appropriate  in  transient  response  prediction  of  s  shock 
wav# excited  infinitely  long  cylindrical  ahell  122!  Thee*  two  values 
of  g  are  appropriate  to  apherical  and  cylindrical  wave  curvature 
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respectively.  tine*  the**  era  the  dominant  acoustic  waveforms  nnr 
ui*  Wet  eurface  of  lhaa*  two  idaauaed  atmctural  bounders** 

Th#  level  of  acoustic  ducmuaujn  iflutd  boundary  element 
densityi  is  an  issus  impacting  the  performance  of  both  DAA2 

form*  In  DAA .  th#  curvature  parameter  l  l«  an  additional 

consideration  that  is  handled  nauonstly  in  DAA'*'  Th*  ability  of 

DAA2  acoustic  fluid  modils  wo  represent  the  fluid  boundary  aigen 
mode*  of  *«ma  aimple  radiating  surface  gt ometrie#  li  of  mureet 
hare  The  above  two  DAA2  versions  can  be  wnttan  in  uncoupled 
mode  form,  these  equations  will  be  used  to  derive  implied  cor- 
responding  specific  acoustic  impedance  functions  for  which  exact 
solutions  exist  for  acoustically  uncoupled  geometrtaa  The  model 
DAA2  equations  era  specialisations  of  equations  (All  end  (A2l. 


DAA7  mnp  ♦  pce^p  ♦  eCui’"'sr.P  “ 

(0(11,11  a  ■»  ”  m,tul.  u.  '*  -  ft  ^ 

DAA  c  (A5i  abova  with  w  :  -  pc  /  -  |  -  c*- 

‘  r' 

substituted  for  w'm 


(A5i 


1A6I 


Geers  has  shown  (16|  that  g  is  related  to  fluid  boundary  atgen 
values  If  1,,  are  the  eigenvalues  of 

m,  ur.  «  i.a^u.  iA?i 

then  the  mode-derived  <*'*'  must  be  restricted  to  the  ran*# 

0  q  <*•  m  <  .  i  Afli 


where 


*r.  “  *n 


These  considerations  given  that  1.  =  mr  a,,,  establish  the 
terms  .  ~  as  shown  in  equation  (ASi 

The  specific  acoustic  impedanoe*  mnlied  by  equations  iA8i 
and  iA6i  (or  radiator  surfaces  demonstrating  uncoupled  modes  can 
!-e  derived  by  apply. ng  [.apiece  transform*  to  (ASi  and  (A6i  Th# 
results  ire 


OAA 


i,  s*  •  pea.  s  *  pcu.  ""  s.  I  p  « 


jpcm.s  *  pcu.  _  m,.sl  u 
OAA  '  iAIOi  above,  with  u  "  replacing  w'J" 


If  velocity  u  is  taken  to  be  modal  velocity  u,  then  equations 
■Ailn  and  ‘Al  !•  can  be  used  to  define  complex-valued  modal 
impedances 


Z.  DAA 


P 

■j, 


pcim.  s  '  e  ,  mrsi 
m.  s*  ♦  «a,.s  *  ecu  '7'  a. 


Since  th*  transform  variabla  s  becomes  >w  for  harmonic  time 
dapendence  of  all  variabla*.  aquations  fAlSi  and  iA13l  become, 
after  torn*  algabra  and  ua*  of  some  new  variables 

"VU,  -  1,. 

°  e  *  *  ~  *rA 


Zn.  DAA  - 


Z_.DAA 


(Cgl^Ju  -  ijw*  ] 


te 


|  'pc i*g  ♦  jcinK  -  !,*«* 


cOr>,.ju  -  l!u* 

L  «c*°n  *  -  **“*  J 


(AlAi 


iA  I5i 


After  further  complex  arithmetic,  and  separation  into  real  and 
imaginary  part*,  the  following  modal  impedance  functions  par 
tlcular  to  tha  two  DAA2  forma  are  obtained 


Re2..  DAA' 


ImZ,..  DAA 
>** 


a*2n.  DAA  1  - 


»c 


i  (E)' 

l  r  n*-./ 1  K/  ~ 


} 


-Alfil 


>A8)  lmZn.  DAA’V 


pc 


«c 


(“■  ‘  r*  (“) )  •  (id) 
r  ©'(?-)’ 


iAITi 


i  A  10) 


(Alii 


where  -  mK  e„.  and  O,.  »  ip  - 

These  axpressions  depend  on  flu.d  bcundary  e'ganvslues  l,.. 
which  are  equivalent  to  modal  accessions  to  ir.orua  per  unit 
radiating  surface  area  mr  *„.  and  the  radiated  wave  curvature 

parameters  g  iDAA  ^'  ,  or  2,  iDAA  '?.*»  It  should  be  pointed  out 

that  fiuic  Boundary  eigenvalue*  correspond  to  low  frequency 
asymptotic  limit*  of  accession  to  inerut  per  unit  anr*  Kxact  ex 
praation*  of  rriodai  aceeasion*  to  insrua  for  the  sphere  and  th*  in 
finite  cylinder  are  given  u-.  Junger  arid  Felt  |£3;  wnr>  »ux-  prest-M 
low  frequency  asymptotic  expressions  Kor  axisy.mmeiric  nu.-les 
the  low  frequency  limn  of  t.  for  a  spherical  surface  :s 


ejr  , 

tr  ispherei  =■  ikn  <<  2n-l 

r‘  n*  l 
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X  DAA  »Ai2>  «bov«  with  wr  ;  r-ipUem^  »«.  | 


i  A  13» 


DTNSRDC  ISSUES  THREE  TYPES  OF  REPORTS 


1.  DTNSRDC  REPORTS,  A  FORMAL  SERIES,  CONTAIN  INFORMATION  OF  PERMANENT  TECH¬ 
NICAL  VALUE.  THEY  CARRY  A  CONSECUTIVE  NUMERICAL  IDENTIFICATION  REGARDLESS  OF 
THEIR  CLASSIFICATION  OR  THE  ORIGINATING  DEPARTMENT. 

2.  DEPARTMENTAL  REPORTS,  A  SEMIFORMAL  SERIES,  CONTAIN  INFORMATION  OF  A  PRELIM¬ 
INARY,  TEMPORARY,  OR  PROPRIETARY  NATURE  OR  OF  LIMITED  INTEREST  OR  SIGNIFICANCE. 
THEY  CARRY  A  DEPARTMENTAL  ALPHANUMERICAL  IDENTIFICATION. 

3.  TECHNICAL  MEMORANDA,  AN  INFORMAL  SERIES,  CONTAIN  TECHNICAL  DOCUMENTATION 
OF  LIMITED  USE  AND  INTEREST.  THEY  ARE  PRIMARILY  WORKING  PAPERS  INTENDED  FOR  IN¬ 
TERNAL  USE.  THEY  CARRY  AN  IDENTIFYING  NUMBER  WHICH  INDICATES  THEIR  TYPE  AND  THE 
NUMERICAL  CODE  OF  THE  ORIGINATING  DEPARTMENT.  ANY  DISTRIBUTION  OUTSIDE  DTNSRDC 
MUST  BE  APPROVED  BY  THE  HEAD  OF  THE  ORIGINATING  DEPARTMENT  ON  A  CASE-BY-CASE 
BASIS. 


